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A is  a general  n x n coefficient  matrix^ 
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INTRODUCTION 

The  general  test  procedure  is  to  apply  each  candidate  routine  to  a 
number  of  test  problems  (with  known  answers)  and  to  record  the  resulting 
accuracy  and  execution  time.  It  is  assumed  that  each  candidate  routine 
can  solve  the  problem 

AX  - B (1) 

where 

A is  a general  n x n coefficient  matrix 
B is  a general  n x m right-hand-side  matrix 
X is  an  n x m solution  matrix 

If  B is  set  equal  to  the  n x n identity  matrix  X (so  that  m ■ n) , then 
X will  be  an  approximation  to  A-1.  The  accuracy  achieved  can  be  measured 
by  computing  the  norm  of  the  residual  matrix 

R - AX  - X (or  R - XA  - X) 

or  better,  by  computing  the  norm  of  the  error  matrix 

E - X - A-1 . 

This  approach  is  used  in  order  to  avoid  having  to  solve  Eq.(l)  repeatedly 
using  various  right-hand  sides  a for  each  test  matrix  A.  Execution  time 
is  measured  by  the  system  clock. 
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a “ “T  11*11  c [4,5] 

n s 

f - -jj-  11*11  F [4,5] 

n>  m ~~  II*  II  „ [5] 

n n 

<J  m "T  ll£Tll  _ [4] 

n S 

Newman  and  Todd  [5]  caution  that  matrices  A and  X exist  for  which 
AX  - I Is  almost  0 while  elements  of  XA  - I are  arbitrarily  large.  In 
their  tests  AX  - I and  XA  - I varied  by  as  much  as  three  orders  of 
magnitude.  Lietzke  et  al.  [4]  report  that 

1.  a < f for  each  test  matrix  and  each  routine  tested. 

2.  a and  f are  not  reliable  estimators  of  q. 

3.  The  following  estimator  for  11*11^  seems  to  work  for  well-conditioned 
matrices: 

IIxrII 

00 

1 “ 1 - 11*11 

00 

where  R - AX  - I. 

Measures  of  error  recommended  here  are  listed  below.  The  Frobenius 
norm  is  recommended,  since  (like  the  Euclidean  matrix  norm)  it  is  com- 
patible with  the  Euclidean  vector  norm  and  Invariant  under  unitary 
transformations. 


Type  of  Error 


Computation 


Actual  relative  error 

Actual  absolute  error 

Estimated  absolute  error 


11*11 , 


ne  II  A"  »ll 


F 


ll**llF 

ne(l  - ll*ll 


F 


Residual  error  11*11  _ 

nt  r 

where  B » X - A-1  and  * ■ AX  - I.  The  above  calculations  must  be  done 
in  double  precision. 
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The  infinity  norm  may  also  be  used — as  a check  on  the  computations, 
and  to  test  the  conjecture  of  Lietzke  et  al.  The  normalization  factor 

~ (where  £ ■ 2 28  s 1.5  x 10-8  is  the  single-precision  machine 
'infinitesimal'  on  UNIVAC  1100  Series  computers)  is  used  to  scale  calcu- 
lated errors  to  unity.  The  normalization  factor  — is  used  to  scale  the 

n 

calculated  errors  to  the  same  units  as  the  elements  of  the  error  matrix. 


TEST  OUTPUT 


For  each  candidate  routine  and  each  test  matrix  the  following  summary 
output  is  needed. 

Heading  identifying  routine 
Name  of  test  matrix 
Order  of  test  matrix 
Value  of  error  return  flag 
Time  for  solution 

Logarithm  of  condition  number  of  test  matrix 

Actual  relative  error 

Actual  absolute  error 

Estimated  absolute  error 

Residual  error 

In  addition,  the  following  detailed  output  should  be  available  upon 
request. 

Test  matrix  A 
Calculated  inverse  x 

Exact  inverse  A-1 

Error  matrix  E - x - A~l 

Residual  matrix  R m AX  - I 


TEST  MATRICES 


The  teat  matrices  recommended  for  use  are  listed  below,  together 
with  a summary  of  their  properties.  These  test  matrices  were  culled 
from  (3,7).  Each  family  of  test  matrices  (except  Wilkinson)  is  defined 
for  any  order  n.  The  test  matrices  and  their  exact  inverses  are  easily 
computed.  The  test  matrices  (except  Newman-Todd)  all  have  integer 
elements  and  span  the  range  from  well-conditioned  to  very  ill-conditioned. 
One  family  (the  Pei  matrices)  have  a parameter  that  allows  one  to  vary 
the  condition  number.  Subroutines  for  generating  the  test  matrices  and 
their  Inverses  are  documented  in  Appendix  B. 
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Matrix 

n 

g 

Condition 

Form 

Wilkinson 

6 

good 

nonsymmetrlc 

inverse  Hilbert*5 

small 

very  bad 

symmetric 

Newman  and  Todd° 

all 

good 

orthogonal,  symmetric 

Rutishauser*5 

small 

very  bad 

lower  triangular 

Peic 

all 

variable 

symmetric  and  positive 
definite 

Givens*5 

all 

fair 

symmetric  and  positive 

definite 


b Precise  condition  numbers  are  given  In  Appendix  A. 

Used  In  study  by  Newman  and  Todd  [5] . 

C Used  by  Newman  and  Todd  15]  for  a ■ n and  by  Lletzke  et  al. 
[4]  for  a - 1. 


Wilkinson  Matrix 


This  matrix  (#3.4  in  [3])  is  for  initial  program  checkout, 
well-conditioned  and  nonsymmetrlc. 


It  is 


” 1 

0 

0 

0 

0 

f 

1 

1 

0 

0 

0 

-1 

W 

-1 

1 

1 

0 

0 

1 

1 

-1 

1 

1 

0 

-1 

-1 

1 

-1 

1 

1 

1 

^1 

-1 

1 

-1 

1 

-1_ 

"~16 

8 

-4 

2 

-1 

1 

-1  1 

0 

16 

8 

-4 

2 

-2 

U m — 

32 

0 

0 

16 

8 

-4 

4 

0 

0 

0 

16 

8 

-8 

0 

0 

0 

0 

16 

16 

Hilbert  Matrix 

_16 

-8 

4 

-2 

1 

-1 

These  matrices  (#3.8 

in 

(31) 

are 

very 

matrix  H of  order  n Is  defined  by 
n J 


*»  " <V 


S 


!* 


•4J 
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where 


(i  , j ■=  1,2 , . . . ,n) 

so  badly  conditioned,  Hilbert  matrices  are  often 
used  to  test  matrix  inversion  routines.  But  this  must  be  done  properly. 

Inverses  of  the  Hilbert  matrices  T ■ H 1 are  used  as  the  test  matrices, 

n n 

since  all  elements  of  T are  integers.  Reference  [3]  lists  the  inverse 

n 

Hilbert  matrices  for  n ■ 2, 3,. ..,10.  For  testing  purposes  it  probably 

suffices  to  try  n = 3,  5,  and  7.  Integer  overflow  will  occur  on  the 
UNIVAC  1110  for  n * 8.  Some  facts  about  Hilbert  matrices,  based  on 

information  in  [2],  are  summarized  below.  Here  C (H  ) denotes  the 

2 n 2 

condition  number  of  (see  Appendix  A)  and  £(7^)  denotes  the  largest 

element  of  r . 

n 


ij  i+j-1 
Because  they  are 


Facts  About  Hilbert  Matrices 


n 

CAH  ) 

2 n 

""A 

Hr  II 

n 2 

*(V 

2 

1.93(+1) 

1.27 

1.52(+1) 

1 . 20(+l) 

3 

5. 24 (+2) 

1.41 

3. 72(+2) 

1.92 (+2) 

4 

1. 55 (+4) 

1.50 

1.03 (+4) 

6. 48 (+3) 

5 

4. 77 (+5) 

1.57 

3. 04 (+5) 

1.79(+5) 

6 

1.50(+7) 

1.62 

9. 24 (+6) 

4. 41 (+6) 

7 

4.75(+8) 

1.66 

2.86(+8) 

1.33 (+8) 

8 

1.53(+10) 

1.70 

9. 00 (+9) 

4. 25 (+9) 

9 

4.93(+ll) 

1.73 

2. 86 (+11) 

1.22 (+11) 

10 

1.60 (+13) 

1.75 

9. 15 (+12) 

3. 48 (+12) 

Newman-Todd  Matrix 

These  matrices  (((3.11  in  [3])  are  orthogonal  and  symmetric,  but  not 
positive  definite.  They  are  defined  by 

' (V’  aa  ' V3rain  iu 

These  matrices  are  well-conditioned,  but  costly  to  generate  (because  of 
calls  to  DSIN  and  DSQRT) . 


5 


I 
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Rutishauser  Matrix 


Theae  matrices  (#3.15  in  [3]  are  ill-conditioned  and  lower  triangular. 
They  are  defined  by 


-2 

-3 

-4 


0 

0 

1 

3 

6 


0 

0 

0 

-1 

-4 


0 

0 

0 

0 

1 


0 

0 

0 

0 


-n+1 


(-Dn_1J 


The  elements  of  R^,  except  for  sign,  are  the  numbers  in  Pascal's  triangle, 


that  is,  binomial  coefficients.  R is  its  own  inverse. 

n 


Pei  Matrix 


The  Pei  matrices  ([3,  p.  18]  and  #18  in  [7])  are  symmetric  and 
positive  definite.  They  are  defined  by 


Va)  “ 


■ij 


a + 1 


i “ 3 


i * 1 


where  a > 0.  Note  that  as  a -*■  0,  P becomes  singular.  Eigenvalues  of  P 
are  X - a + n (simple)  and  X - a (multiplicity  n - 1) . Thus,  for  a as  n. 


P is  well-conditioned;  but  for  a a 0,  P is  very  ill-conditioned. 
n n J 


Recommended  choices  of  the  Pei  parameter  a are 


a - 64e,  1,  n 


where  E ■ 2 26  a 1.5  x 10  8, 


To  compute  the  inverse  of  P (a) , let  J be  the  matrix  with  all  elements 

n 


equal  to  1.  Then  (dropping  the  subscript  n) 
P - a!  + J. 


rik" jar-' 


— : — = — 
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2 “1 

Since  J ■ nJ,  we  expect  P = Q to  be  linear  in  J, 

0 m bl  + cJ. 

Computing  the  product  PQ  we  find 

1 - PQ  - (aJ  + J)(bl  + cJ)  - abl  + (jb  + ac  + nc)J. 


Equating  coefficients, 


ab  - 1 


a 


b + ac  + nc  m 0 . 


a+n  a ( a+n ) 


Hence, 


c ■ i (J  - ss J) 


p"1  • (ijj). 


a + n - 1 
a (a+n) 


1 

a (a+n) 


i - j 


i =£  j . 


Givens  Matrix 

The  Givens  matrices  (#8  in  (7J)  are  poorly-conditioned,  symmetric, 
and  positive  definite.  They  are  defined  by 

“l  1 1 ...  1 

1 3 3 ...  3 

1 3 5 ...  5 


• • • • 


1 3 5 ...  2n-l 


Thus,  Gn  - i9ij)  where  gr^.  - 2 min(i,j)  - 1.  The  inverse  is  given  by 


I 

I 
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Note  that  the  inverse  Givens  matrix  is  tridiagonal. 


1 


J 


/<£ 
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Appendix  A. 


MATRIX  NORMS  AND  CONDITION  NUMBERS 


The  more  common  matrix  and  vector  norms  are  defined  here,  without 
much  explanation.  For  a more  complete  discussion  of  matrix  and  vector 
norms  and  their  properties,  see  [1,6].  Vector  norms  in  common  use  are 


Ixll.  - 2. |*. | 

1 j. ' i 1 


Ixll 2 - ^EjX* 


(Euclidean) 


llxll^  - maxj  |xi| 


(Chebyshev) 


These  are  particular  instances  of  the  p-norm 


Ixll  - EJ* 


Any  vector  norm  induces  a matrix  norm  according  to  the  following  relation: 


IIaII  ■ max..  ..  . lUxll 
llxll-1 


(A.l) 


where  the  maximum  is  taken  over  all  x such  that  llxll  ■ 1.  Matrix  norms 
induced  in  this  way  are  compatible  with  the  underlying  vector  norm,  in 
the  sense  that 


llAxll  < IUII  • llxll 


(for  all  x) . 


The  matrix  norms  Induced  by  the  vector  p-norms  (p  - 1,  2,  ®)  are 


IUII  J - many 


(column  sum) 


IUII  2 " ^largest  eigenvalue  of  ATA 


(Euclidean) 


(row  sum) 


/S  H 

r— — — * ? 
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Matrix  norms,  besides  those  induced  by  vector  p-norms,  are  also  used. 
Among  the  most  useful  are 


I'Al's 

H AH  p 


E.S.I a.  ,| 

1 j'  ij 1 


n max . max  . 
i J 


(sum  norm) 
(Frobenius  norm) 
(max  norm) . 


It  is  easily  proven  that  IUII  is  compatible  with  llxll  , that  IUII  is 

S 1 r 

compatible  with  llxll 2,  and  that  lUII^  is  compatible  with  llxll oo.  Also, 

i II  All  c <IUII  < IUII  < nil  All 

n S 00  M S 

7 Ha  II  <llAll,  < IUII  < nil  All 

n o 1 M S 

so  that  the  choice  of  a norm  is  not  too  critical  in  finite-dimensional 
linear  spaces  (that  is,  for  matrices). 

The  condition  number  of  a matrix  for  a given  matrix  norm  is  defined 
by 

c(a)  - IUII  • IU-1  II.  (A. 2) 

From  this  definition  it  is  clear  that  c(A  1 ) * c(A) . The  condition 
number  recommended  in  this  report  is  derived  from  the  Frobenius  matrix 
norm 

cf(a)  - I|aIIf  • IU‘1IIF. 

The  Frobenius  condition  number  should  be  a good  approximation  to  the 
Euclidean  condition  number 

c2(a)  “ ^largest  eigenvalue  of  A^X/smallest  eigenvalue  of  ATA 


because  of  the  inequalities 

IU II 2 < IIaIIf  < /nil All 2 . 

A specialization  of  the  Euclidean  condition  number,  introduced  by 
von  Neumann  and  Goldstine  in  the  early  days  of  computing,  is  frequently 
found  in  the  literature  (for  example,  (3)).  This  condition  number  Is 
defined  by 
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c (A)  _ I largest  eigenvalue  of  A I 
A I smallest  eigenvalue  of  A I 

It  is  sometimes  called  the  'spectral'  condition  number,  since  the 
spectral  radius  of  A is  the  magnitude  of  the  largest  eigenvalue  of  A. 
When  A is  symmetric,  the  Euclidean  and  spectral  condition  numbers  agree 

c2(A)  - c^(A). 

Spectral  condition  numbers  for  some  of  the  test  matrices  discussed 
in  this  report  are  given  below.  These  are  quoted  from  [3J . 

matrix  spectral  condition  number 

Newman-Todd 

Rutishauser 

Pei 

Givens 

Condition  numbers  computed  from  Eq.  (A. 2)  for  various  matrix  norms  will, 
because  of  the  inequalities  relating  the  matrix  norms,  all  be  of  the 
same  order  of  magnitude  as  the  spectral  condition  numbers  quoted  above. 

For  any  condition  number  c(A) , Forsythe  and  Moler  [2]  argue  that 

log j 0 c(A) 

is  approximately  the  number  of  significant  digits  lost  when  calculating 

A . A little  arithmetic,  based  on  the  condition  numbers  quoted  above, 
yields  the  following: 


1 

4 n £.n  2 


1 + - 

a 

- {—) 
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matrix 

order 

lo«10  c 

Inverse  Hilbert 

3 

2.72 

5 

5.68 

7 

7.18 

Newman-Todd 

all 

0 

Rutishauser 

3 

3.61 

5 

6.02 

7 

8.43 

Pei  a - n 

all 

0 

a - 1 

10 

1. 

50 

1.70 

100 

2. 

a “ 64e 

10 

7.02 

50 

7.72 

100 

8.02 

a - e 

10 

8.83 

50 

9.53 

100 

9.83 

Givens 

10 

2.71 

50 

4.10 

100 

4.71 

From  the  above  table  it  is  apparent  that  condition  numbers  for 
families  of  dense  matrices  of  orders  up  to  several  hundred  can  be 
classified  as  follows: 

Condition  Functional  Form  Example 

Good  constant  Nevman-Todd,  Pei  (a  ■ n) 

Fair  linear  in  n Pei  (a  =*  1) 

quadratic  in  n Givens 

Very  bad  exponential  in  n Rutishauser,  Hilbert 

-20 

Note  too  that  Pei  (a  - 64e  *2  ) presents  a more  interesting  test 

problem  than  Pei  (a  ■ e) . 
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Appendix  B 

PROGRAM  PACKAGE  FOR  TESTING  LINEAR  EQUATION  ROUTINES 

A program  package  for  testing  linear  equation  routines  was  written 
according  to  the  specifications  of  this  report.  It  consists  of  the 
following: 

Main  Program:  LEQTST  - overall  logic 

Subprograms:  LTxxxx  - generates  test  matrix  and  its  exact 

inverse 

LTEST  - calls  candidate  routine  to  calculate 
inverse  and  measures  execution  time 

LTERR  - computes  error  matrix,  residual  matrix 
and  four  measures  of  error 

FNORM  - returns  Frobenius  norm  of  a matrix 

LTPRNT  - prints  summary  test  results 

LTRITE  - writes  all  matrices,  if  detailed  printout 
requested 

MXRITE  - writes  a single  matrix 

Program  documentation  for  these  program  elements  is  contained  in  this 
Appendix,  and  conforms  to  the  documentation  standards  of  [8], 


I 
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FNORM-1 


NAME 

PURPOSE 


USAGE 


ACCESS 

ERRORS 

REMARKS 

PROGRAM  INFO 


TESTING 

REFERENCE 


FNORM 

To  return  Frobenlus  norm  of  a matrix 


F = FNORM  (A,  N,  NN) 

A The  matrix  - ( N , N ) 

N Order  of  matrix  A 

NN....Row  dimension  of  A 


Input 

Input 

Input 


LIB  NWC*MATHLIB. 

None 

Although  A is  single  precision,  the  computation  of  FNORM 
is  done  in  double  precision. 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 

Filename 

Element/Vers 

Storage 

Timing 

Consultant 


UNIVAC  1110 
FORTRAN  V 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 

Certified,  Fully  Supported  by  CCF 
FNORM 

DSQRT,  NERR3$ 

NWC*MATHLIB 
FNORM 
74  words 
Unknown 

L.  W.  Lucas,  Ext.  3561 


Output  from  FNORM  was  checked  against  values  computed  by 
hand. 

NWC  TM  2410 


15 


TM  2410 


LEQTST-1 


NAME 

PURPOSE 

USAGE 


ACCESS 

REMARKS 


LEQTST 

To  test  a linear  equation  routine  by  using  it  to  compute 
inverses  X of  tests  matrices  A having  known  inverses. 

The  resulting  error  E = X - A-1  and  residual  R m AX  - I 
are  used  to  compute  several  measures  of  computational 
error. 

Main  program.  There  are  no  inputs  to  LEQTST.  User  options 
involve  the  setting  of  three  program  variables. 

NN Parameter  variable  which  controls  the 

maximum  size  of  test  matrices  used. 

DETAIL. .. .Logical  variable  which  if  .TRUE,  causes  all 
matrices  involved  - A,  X,  A-1,  E,  R - to  be 
printed  out,  as  well  as  summary  test  informa- 
tion. 

IN Integer  variable  which  specifies  the  logical 

unit  number  for  output. 

LEQTST  is  set  up  to  use  the  following  test  matrices. 


matrix 

order 

Wilkinson 

6 

Inverse  Hilbert 

3, 

5, 

7 

Newman-Todd 

5, 

10, 

50, 

100 

Rutishauser 

5, 

10, 

15, 

20 

Pei,  a = 64e 

5, 

10, 

50, 

100 

a » 1 

5, 

10, 

50, 

100 

a = n 

5, 

10, 

50, 

100 

Givens 

5, 

10, 

50, 

100 

The  maximum  order  actually  used  is  controlled  by  the 
value  of  NN. 

IN  4033519*NM-BENCH. LEQTST 

The  size  of  test  matrices  used  in  LEQTST  is  limited  only 
by  the  available  core  storage.  The  design  of  LEQTST  does 
not  provide  for  buffering  to  mass  storage.  LEQTST  is  not 
intended  to  test  sparse  matrix  routines.  The  test  matrices 
used  are  believed  to  be  representative  of  moderately-sized 
dense  matrices. 


TM  2410 


LEQTST-2 


r1 

i 


PROGRAM  INFO 


TESTING 

REFERENCE 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 


Filename 

Element/Vers 

Storage 

Timing 

Consultant 


UNIVAC  1110 
FORTRAN  V 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 
Unsupported 
None 

LTHDG,  LTWILK,  LTEST,  LTERR,  FNORM, 
LTPRNT,  LTRITE , LTHILB,  LTNEWT , LTRUTH, 
LTPEI,  LTGIVN,  NINTR$ , ALOG,  NST0P$ 
4033519*NM-BENCH • 

LEQTST 

51080  words  (AW  ■ 100) 

103  seconds  (NN  - 100) 

None 


The  component  subprograms  for  LEQTST  were  all  handchecked. 
NWC  TM  2410 


TM  2410  LTERR-1 

NAME  LTERR 

PURPOSE  To  compute  error  and  residual  matrices  and  four  measures 


of  computational  error. 

USAGE  CALL  LTERR(A,  X,  AI,  E,  R,  N,  NN,  ERR) 

A. ...Test  matrix  - (W,  N ) Input 

X. .. .Computed  inverse  of  A Input 

AI... Exact  inverse  of  A Input 

E.... Error  matrix,  E « X - AI  Output 

R. .. .Residual  matrix,  R - AX  - I Output 

N.... Order  of  matrices  A , X,  AI,  E,  R Input 
NN . . . Row  dimension  of  A , X,  AI , E,  R Input 
ERR.. Array  of  error  measures  output 

ACCESS  IN  4033519*NM-BENCH. LTERR 

ERRORS  None 


REMARKS  Although  the  arrays  A,  X,  AI , E and  R are  all  single 

precision,  the  computations  within  LTERR  are  performed  in 
double  precision.  The  measures  of  error  returned  in  ERR 
are  defined  in  the  reference. 

PROGRAM  INFO  ' 

Machine  UNIVAC  1110 

Language  FORTRAN 

Author  L.  W.  Lucas,  Code  4033  NWC 

Date  14  March  1975 

Status  Unsupported 

Entry  Names  LTERR 

External  Refs  FNORM,  DSQRT , NERR3$ 

Filename  40335 19*NM-BENCH.  , 

Element /Vers  LTERR 

Storage  333  words 

Timing  unknown 

Consultant  None 

TESTING  Output  from  LTERR  was  checked  against  values  computed  by 

hand. 

REFERENCE  NWC  TM  2410 


TM  2410 


LTEST-1 


I 


I 

I 

( 


NAME  LTEST 

PURPOSE  To  call  candidate  routine  for  calculating  Inverse  of  test 

matrix  and  to  measure  execution  time 

USAGE  CALL  LTEST  (A,  X,  N,  NN,  TIME,  WK,  C,  IER) 


A • • • • 

..Test  matrix  - (N,  N) 

Input 

X. . . . 

..Computed  inverse  of  A 

Output 

N 

..Order  of  matrices  A,  x 

Input 

NN... 

..Row  dimension  of  A,  X 

Input 

TIME. 

..Time  to  compute  inverse 

Output 

WK... 

..Work  space  of  length  N*N 

Scratch 

C 

..Work  space  - (N,  N) 

Scratch 

IER. . 

..Error  return  flag 

Output 

ACCESS  IN  40335 19*NM-BENCH. LTEST 

ERRORS  No  errors  may  occur  within  LTEST  itself.  IER  returns  the 

value  of  an  error  flag  from  the  candidate  routine  (if  any). 

REMARKS  The  input  matrix  A is  copied  into  the  work  space  C,  which 

is  passed  to  the  candidate  routine,  just  in  case  the  rou- 
tine destroys  any  coefficient  matrix  given  to  it.  LTEST 
will  obviously  have  to  be  rewritten  for  each  candidate 
routine  tested.  SETCLK  and  LKCLKS  are  entry  points  to  the 
system  clock  routine.  LEQT1F  is  the  routine  being  tested. 

PROGRAM  INFO 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 

Filename 

Element/Vers 

Storage 

Timing 

Consultant 

REFERENCE  NWC  TM  2410 


UNIVAC  1110 
FORTRAN 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 
Unsupported 
LTEST 

SETCLK,  LEQTIF,  LKCLK,  NERR3$ 
4033519*NM-BENCH . 

LTEST 
167  words 
unknown 
None 


21 
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1 


f 

r 

L 

TM  2410  LTGIVN-1 

i 

NAME  LTGIVN 

PURPOSE  To  generate  the  Givens  test  matrix  of  order  N and  its 

exact  inverse. 

Outpu 
Output 
Input 
Input 


USAGE  CALL  LTGIVN  (G,  GI,  N,  NN) 

G Test  matrix  - (N,  N ) 

GI . . . . Exact  inverse  of  G 

N Order  of  G 

NN....Row  dimension  of  G 

ACCESS  LIB  NWC*MATHLIB. 

ERRORS  None 


? 


it 


REMARKS  None 

PROGRAM  INFO 


TESTING 

METHOD 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 

Filename 

Element/Vers 

Storage 

Timing 

Consultant 


UNIVAC  1110 
FORTRAN 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 

Certified,  Fully  Supported  by  CCF 

LTGIVN 

NERR3$ 

NWC*MATHLIB 
LTGIVN 
198  words 
unknown 

L.  W.  Lucas,  Ext.  3561 


Output  from  LTGIVN  was  handchecked  against  Gregory  and 
Karney. 

The  Givens  test  matrix  and  its  inverse  are  defined  by 

“l  1 1 ...  1 ~ 

1 3 

G - 

n 

1 3 5 ...  2n-l 


I 

t 


i 


r 


3 -1 
-1  2 


2 

1 


1 

1 


| 


REFERENCE  Gregory  and  Karney.  A Collection  of  Matrices  for  Tenting 
Computational  Algorithms.  Wiley-Interacience , 1969. 


| f 

TM  2410  LTHILB-2 

REFERENCES  Gregory  and  Karney 

A Collection  of  Matrices  for  Testing  Computational 
Algor i thms 

Wiley-Interscience,  1969 
Forsythe  and  Moler 

Computer  Solution  of  Linear  Algebraic  Systems 
Prentice-Hall,  1967 


TM  2410 


LTNEWT-1 


NAME 

PURPOSE 

USAGE 

ACCESS 

ERRORS 

REMARKS 

PROGRAM  INFO 


TESTING 

METHOD 

REFERENCE 


LTNEWT 

To  generate  Newman-Todd  test  matrix  of  order  N and  its 
exact  inverse. 

CALL  LTNEWT  (A,  AI,  N,  NN) 


A Newman-Todd  test  matrix  - (N,  N)  Output 

AI.... Exact  inverse  of  A Output 

N Order  of  A Input 

NN . . . . Row  dimension  of  A,  AI  Input 


LIB  NWC*MATHLIB 
None 

Since  A is  orthogonal,  AI  is  set  equal  to  the  transpose 
of  A.  The  computations  are  done  in  double  precision, 
since  the  elements  of  A are  not  exactly  representable  in 
single  precision. 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 

Filename 

Element/Vers 

Storage 

Timing 

Consultant 


UNIVAC  1110 
FORTRAN 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 

Certified,  Fully  Supported  by  CCF 
LTNEWT 

DSQRT , DSIN,  NERR3$ 

NWC*MATHLIB 
LTNEWT 
142  words 
unknown 

L.  W.  Lucas,  Ext.  3561. 


Output  from  LTNEWT  was  handchecked  against  Gregory  and 
Karney. 

The  Newman-Todd  test  matrix  A is  defined  by 


A - (a^). 


Lj 


Gregory  and  Karney 

A Collection  of  Matrices  for  Testing  Computational 
Algorithms 

Wiley-Interscience,  1969 


TM  2410 


LTPEI-1 


NAME  LTPEI 

PURPOSE  To  generate  Pel  test  matrix  of  order  N with  parameter  S 

and  its  exact  inverse. 


USAGE  CALL  LTPEI  (P,  Q,  N,  NN,  S) 


p... 

..Pei  test  matrix  - (N,  N) 

Output 

Q... 

. . Exact  inverse  of  P 

Output 

N... 

. .Order  of  P , Q 

Input 

NN.. 

..Row  dimension  of  P,  Q 

Input 

S... 

..Pei  parameter 

Input 

ACCESS  LIB  NWC*MATHLIB. 

ERRORS  None 

REMARKS  The  elements  of  P are  integers  and  are  exact  to  single 

precision.  The  elements  of  Q are  fractions  and  are  com- 
puted in  double  precision.  The  parameter  S can  be  used 
to  vary  the  condition  number  of  P. 

PROGRAM  INFO 


TESTING 

METHOD 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 

Filename 

Element/Vers 

Storage 

Timing 

Consultant 


UN1VAC  1110 
FORTRAN 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 

Certified,  Fully  Supported  by  CCF 

LTPEI 

NERR3$ 

NWC*MATHLIB 
LTPEI 
183  words 
unknown 

L.  W.  Lucas,  Ext.  3561 


Output  from  LTPEI  was  handchecked  against  Gregory  and 
Karney . 

The  formulas  defining  P and  0,  taken  from  Gregory-Karney , 
are 


p ■ <V 


"a 


1 + s 
1 


i - j 
i + J 


0 “ (*ij> 


qH 


s + n - 1 
s(s  + n) 

_ 1 
s(s  + n) 


where  s is  the  Pei  parameter. 


i - J 
i * Jf 


29 


1 

I 


mm  .irUiifflffr""1 


TM  2410 


LTPEI-2 


REFERENCE 


Gregory  and  Karney 

A Collection  of  Matrices  for  Testing  Computational 
Algor i thms 

Wiley-Inter8cience,  1969 


30 


.. , •. 


I 


I 


I 


T 


TM  2410 


LTPRNT-1 


NAME  LTPRNT 

PURPOSE  To  print  summary  test  results  for  LEQTST 

USAGE  CALL  LTHDG  (IW) 

IW Unit  number  for  printout 

CALL  LTPRNT  (IW,  NAME,  N,  IER,  TIME,  COND,  ERR) 

Input 
Input 
Input 
Input 
Input 
Input 
Input 


IW Unit  number  for  printout 

NAME... Name  of  test  matrix 

N Order  of  test  matrix 

IER.... Error  return  from  candidate  routine 
TIME... Time  to  compute  inverse 
COND. . .Logarithm  of  condition  number 
ERR. . . .Array  of  error  measures 


ACCESS  IN  4033519*NM-BENCH. LTPRNT 

ERRORS  None 

REMARKS  NAME  is  a 12-character  hollerith  string  giving  the  name 

of  the  test  matrix.  Entry  LTPRNT  is  used  to  print 
summary  test  results.  Entry  LTHDG  is  used  to  print 
headings  for  the  summary  test  results. 

PROGRAM  INFO 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 

Filename 

Element/Vers 

Storage 

Timing 

Consultant 


UNIVAC  1110 
FORTRAN 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 
Unsupported  • ' 

LTHDG,  LTPRNT 

NWDU$ , NI01$ , NI02$,  NERR3$ 

4033519*NM-BENCH 

LTPRNT 

108  words 

unknown 

None 


REFERENCE  NWC  TM  2410 


I 

i 


i 


31 


TM  2410 


LTRITE-1 


NAME  LTRITE 

PURPOSE  To  write  matrices  for  LEQTST 

USAGE  CALL  LTRITE  (IW,  A,  X,  AI,  E,  R,  N,  NN) 


IW.. 

...Logical  unit  number  for  output 

Input 

A. . . 

. .Test  matrix  - (N,  N) 

Input 

X... 

..Computed  Inverse 

Input 

AI.. 

..Exact  inverse 

Input 

E • • • 

..Error  matrix 

Input 

R.  • • 

..Residual  matrix 

Input 

N... 

..Order  of  matrices 

Input 

NN.. 

..Row  dimension  of  arrays 

Input 

ACCESS  IN  4033519*NM-BENCH. LTRITE 

ERRORS  None 

REMARKS  None 

PROGRAM  INFO 

Machine 
Language 
Author 
Date 
Status 
Entry  Names 
External  Refs 
Filename 
Element /Vers 
Storage 
Timing 
Consultant 


UNIVAC  1110 
FORTRAN 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 
Unsupported 
LTRITE 

MXRITE,  NERR3$ 

4033519*NM-BENCH 

LTRITE 

96  words 

unknown 

None 


REFERENCE  NWC  TM  2410 


33 


TM  2410 


LTRUTH-1 


NAME 

PURPOSE 

USAGE 


ACCESS 

ERRORS 

REMARKS 

PROGRAM 


TESTING 

METHOD 


LTRUTH 


To  generate  Rutlshauser  test  matrix  of  order  N and  its 
exact  Inverse. 


CALL  LTRUTH  (R,  RI,  N,  NN,  JR) 

R Rutlshauser  test  matrix  - (W,  N) 

RI.... Exact  Inverse  of  R 

N Order  of  R,  RI 

NN....Row  dimension  of  R,  RI 
JR.... Work  space  array  of  length  N 


Output 

Output 

Input 

Input 

Scratch 


LIB  NWC*MATHLIB. 


None 

The  scratch  array  JR  is  used  to  store  the  jth  column  of  R, 
while  generating  the  j+l-th  column  via  a recursion 
formula.  R is  its  own  inverse,  so  RI  is  just  a copy  of  R. 

INFO 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 

Filename 

Element  Vers 

Storage 

Timing 

Consultant 


UNIVAC  1110 
FORTRAN 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 

Certified,  Fully  Supported  by  CCF 

LTRUTH 

NERR3$ 

NWC*MATHLIB 
LTRUTH 
197  words 
unknown 

L.  W.  Lucas,  Ext.  3561 


Output  from  LTRUTH  was  handchecked  against  Gregory  and 
Karney . 


The  Rutlshauser  test  matrix  is  defined  by 


0 


35 


TM  2410 


LTRUTH-2 


REFERENCE 


The  columns  of  R are,  except  for  sign,  the  diagonals  In 
Pascal's  triangle.  Thus,  the  following  recursion  formula 
can  be  used  to  generate  the  elements  of  R. 


'll'* 


ri-l,  j " ri-l,  j- 1 


i < j 
i > j 


Gregory  and  Karney 

A Collection  of  Matrices  for  Testing  Computational 
Algorithms 

Wiley-Interscience , 1969 


/ 
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TM  2410 


LTWILK-1 


NAME  LTWILK 

PURPOSE  To  generate  Wilkinson  test  matrix  of  order  6 and  its  exact 

inverse. 

USAGE  CALL  LTWILK  (W,  V,  N,  NN) 


W Wilkinson  test  matrix  Output 

V Exact  inverse  of  W Output 

N Order  of  W,  V Output 


NN....Row  dimension  of  W,  V Input 

ACCESS  LIB  NWC*MATHLIB. 

ERRORS  None 

REMARKS  N is  set  equal  to  6 by  LTWILK. 

PROGRAM  INFO 

Machine  UNIVAC  1110 

Language  FORTRAN 

Author  L.  W.  Lucas,  Code  4033  NWC 

Date  14  March  1975 

Status  Certified,  Fully  Supported  by  CCF 

Entry  Names  LTWILK 

External  Refs  NERR3$  1 

Filename  NWC*MATHLIB 

Element /Vers  LTWILK 

Storage  165  words 

Timing  unknown  , 

Consultant  L.  W.  Lucas,  Ext.  3561 

TESTING  Output  from  LTWILK  was  handchecked  against  Gregory  and 

Karney. 

METHOD  The  Internal  arrays  X and  Y are  initialized  to  the  Wilkinson  I 

test  matrix,  and  its  exact  inverse,  respectively,  which  are 
given  in  the  reference.  These  are  then  copied  into  W and 
V,  respectively,  upon  entry  to  LTWILK. 

REFERENCE  Gregory  and  Karney 

A Collection  of  Matrices  for  Testing  Computational 
Algorithms 

Wiley-Interscience,  1969 


TM  2410 


MXRITE-1 


NAME 

PURPOSE 

USAGE 


ACCESS 

ERRORS 

REMARKS 


PROGRAM  INFO 


REFERENCE 


MXRITE 


To  write  out  a general  matrix,  preceded  by  a title. 
CALL  MXRITE  (IW,  A,  N,  M,  NN,  TITLE) 


IW Logical  unit  number  for  printout 

A The  matrix  (N,  M) 

N Rows  in  A 

M Columns  in  A 

NN Row  dimension  of  A 

TITLE.. Title  printed  above  A 

NWC*MATHLIB. 

None 

TITLE  is  an  18-character  hollerith  string  printed  above  A. 
The  matrix  is  printed  out  8 columns  across  the  page. 
Remaining  columns  are  printed  below,  again  8 columns 
across  the  page. 


Input 

Input 

Input 

Input 

Input 

Input 


Machine 

Language 

Author 

Date 

Status 

Entry  Names 

External  Refs 

Filename 

Element/Vers 

Storage 

Timing 

Consultant 

None 


UNIVAC  1110 
FORTRAN 

L.  W.  Lucas,  Code  4033  NWC 
14  March  1975 

Certified,  Fully  Supported  by  CCF 
MXRITE 

NWDU$,  NI01$,  NI02$ , NERR3$ 

NWC*MATHLIB 

MXRITE 

167  words 

unknown 

L.  W.  Lucas,  Ext.  3561 
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